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ABSTRACT 

We present simple analytical expressions for the predicted spectral and temporal 
behavior of the early afterglow radiation from gamma-ray bursts in radiative regimes 
intermediate between the adiabatic and the fully radiative solutions of the blastwave 
hydrodynamic equations. Our expressions are valid as long as the relativistic electrons 
responsible for the observed synchrotron emission are in the fast cooling regime and 
the blast wave is relativistic. We show that even a slight deviation from a perfectly 
adiabatic evolution results in significant changes of the temporal characteristics of the 
afterglow emission. 
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Introduction 



The relativistic blast wave model ( |Meszaros fc Rees 1992| , Rees & Meszaros 1992) has met 
with considerable success in explaining the X-ray, optical and radio afterglows of gamma-ray 
bursts (GRBs). In this model, the impulsive release of ~ lO^'^ - 10^^ ergs of energy (e. g., Kumar 



1999P into a small volume results in the formation of a relativistic blast wave. The initial Lorentz 
factor Tq of material behind the shock wave reaches values of Tq ~ 



100 - 1000 (iBlandford fc 



McKee 1976 , Cavallo &: Rees 1978 , Shemi &: Piran 1990| ). The dominant radiation mechanism 



responsible for the radio - X-ray afterglows from these relativistic blast waves is believed to be 
optically thin synchrotron emission from relativistic electrons accelerated behind the shock front 
QMeszaros fc Rees 1993| , [Katz 1994| , p^avani 1996| ). 
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With rapid follow-up observations of the X-ray, optical, and radio afterglows of several GRBs 
(e. g., [Costa et al.T997i |Piro et al. l"998| , [Metzger et al. 199% [van Paradijs et al. 1997 , Djorgovski[ 
[et al. 1997 , Kulkarni et al. 1998, Frail et al. 1997) it has become possible to determine in great 
detail the spectral and temporal characteristics of these afterglows, which can then be compared 
to the predictions of the relativistic blast wave model ( Palama et al. 1998b[ , [Iwamoto 19"99[ ). 



The temporal evolution of flux levels and characteristic frequencies of the synchrotron 
spectrum from electrons in a relativistic blast wave can be formulated in a very elegant way, 
directly comparable to observations, if the blast wave is either perfectly adiabatic or fully radiative 
( Sari et al. 1998| ). However, while the comparison to afterglow data of several GRBs suggests 
that at least in the late afterglow phase GRB blast waves are well described by the adiabatic 
solution (e. g., Galama et al. 1998b| ), they can obviously not be perfectly adiabatic in order to 
be observable, and they are generally believed to be strongly non-adiabatic in the early afterglow 
phase. 

Analytic estimates similar to the ones found by Sari et al. (1998) have been developed for 
intermediate radiative regimes under the a priori assumption that the blast wave evolution follows 



a self-similar behavior T{r) oc r ^ (Rees & Meszaros 1998, Dermer et al. 1999, Meszaros & Rees 
1999). However, the relation between the deceleration index g and the blast wave energetics and 



radiative properties remains uncertain in these representations because the blastwave dynamical 
equations are not solved self-consistently. 

Self-consistent solutions of the blast wave dynamics under realistic assumptions on the energy 
transfer from protons to electrons and on the magnetic field evolution, calculating the synchrotron 
(and synchrotron self Compton) spectrum, solving for the evolution of the electron distribution 
behind the shock front and the blast wave dynamics simultaneously, have up to now only been 
possible numerically (Chiang & Dermer 1999, Huang et al. 199"9| , Moderski et al. 1999| ). 



In this paper, we show that under certain conditions, if the electrons are in the fast-cooling 
regime and the blast wave is relativistic, the blast wave kinetic equation can be solved analytically, 
yielding a simple, self-consistent analytical representation for the synchrotron emission from 
relativistic blast waves in general radiative regimes. In Section 2, we quote the blast wave kinetic 
equation and present its solution under the above assumptions. The spectral and temporal 
characteristics of the synchrotron radiation from relativistic blast waves are derived in Section 3. 
In Section 4, we discuss implications of a non-zero radiative efficiency for the interpretation of 
observed GRB afterglow characteristics. We summarize in Section 5. In this paper, we neglect the 
effects of synchrotron self-absorption and synchrotron-self Compton scattering. The importance 
of these effects is investigated in detail in a separate paper (Dermer et al., in preparation). 
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2. The kinetic equation for general radiative regimes 

In the relativistic blast wave model it is assumed that a total energy E = 10^^ £'52 erg is 
deposited into a small volume, giving rise to a relativistic blast wave expanding with a bulk 
Lorentz factor Fq. The mass of the initial baryon-loaded ejecta is given by Mq = Eo/iJ'QCP') in 
the comoving frame of the shocked material. As it expands into an external medium of density 
'^extif) oc r~^^ it sweeps up matter at a rate dm = Q.r'^ nf,xt{r)mpdr, where is the solid angle 
element into which the blast wave is expanding. The kinetic equation governing the blast wave 



evolution is then given by (Blandford & McKee 1976, Chiang & Dermer 1999) 



- = -^ (1) 
dm M ' ^ ' 

where M is the total, comoving mass (rest mass plus internal kinetic energy) of the ejecta + 

swept-up material. We assume that a fraction Ce of the swept-up energy per unit time is transferred 

to relativistic electrons behind the shock front, and a fraction Crad of this energy is then radiated 

as synchrotron radiation. Thus, a fraction e = ee^rad of the swept-up energy will be transformed 

into radiation. If the relativistic electrons are in the fast cooling regime (i. e. the synchrotron 

cooling time scale for all electrons is shorter than the dynamical time scale), then the radiative 

efficiency e^ad ~ li and e f« eg may be regarded as constant over the fraction of the blast wave 

evolution during which this condition is met. In this regime, Eq. || can be solved analytically. We 

note that up to now it is not certain whether the condition ~ const, is actually met in realistic 

blastwave scenarios. Investigation of this question would require the detailed treatment of the 

process of energy transfer from protons to electrons (see, e. g., Pohl & Schlickeiser 1999), which is 

beyond the scope of this paper. Here, we assume that eg ~ const, and note that our results may 

not be used if eg or e^ad change significantly during the early afterglow phase. 

As the blast wave propagates through the surrounding medium, its mass (in the comoving 
frame) increases at a rate 

dM ={e + T[l-€]) dm. (2) 

Inserting Eq. ^ into Eq. |l], we find 

-'2 1 \ 2 



M(r) = M„ \—A -^1 . (3) 



. 

Using this relation in Eq. |l|, we have the general solution 



m[r) f (7 + 1 



MoAo 7 (72-1)3/2' 
where 
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and 



The integral on the right-hand-side of Eq. ^ can be solved analytically in the extreme cases of an 
adiabatic (e = 0) and a fully radiative (e = 1) blast wave, yielding 

m{r) _ Fad To 



and 

m{r) ^ Trad + 1 _ Tq -H 

in the adiabatic and radiative limits, respectively. It is straightforward to show that in the 
relativistic limit (Fq > F » 1) and in the deceleration phase (m(r) » Mq/Fq), these solutions 
approach the asymptotic behaviors F^^ oc r~^^~^^/'^ and Frarf(r) cx r~^^~^\ if the external density 
is described by a radial profile 



riextir) = l^—j . (9) 
For general radiative regimes, the integral in Eq. |^ can be solved in the relativistic limit, yielding 

Mo 



, r , m(r) Fn\ '--^ 
Trei{r) = To l + [2-e] \' " ) , (10) 



which has the asymptotic limits 



if (2 - e) m(r) Fq « Mq 
r.e/(r) ^{ ^ ^ (11) 
a if (2 - e) m(r) Fq » Mq 

where 

[3-^]^oro- .2) 
I [2 — e] n,.(, TUpC^ Tq / 

Being interested in the spectral characteristics of the afterglow radiation, we will use the late-time 
asymptote of Eq. |ll| in the remainder of this paper. Integrating over time t in the observer's 
frame, using dt = dr/iT'^c), we find 



a? ct 



and 



where 



m = \-^ (13) 



nt) = a^ {%]\ (14) 



P = (15) 
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and 

„-3 

8-27?-e 2 ^ ^ 

In the limits e = and e = 1, Eqs. [l^ - |l^ reproduce the well-known scaling laws for adiabatic 
and the radiative blast waves, respectively. 



3. Temporal evolution of synchrotron spectra 

For the purpose of illustration, we will derive the relevant scaling laws of the early afterglow 
emission here for the special case t/ = 0, i. e. for a homogeneous external medium of density 
no cm~^. A fraction eg (~ e in the fast-cooling regime) of the swept-up energy is assumed to be 
transferred to a relativistic, nonthermal electron distribution injected into the post-shock region 
with a power-law spectrum of index p {n{j) oc 7~^), where we generally assume 2 < p < 3. 
With the energy normalization determined by the electron acceleration efficiency Ce, the electron 
injection distribution must have a low-energy cutoff at Lorentz factor (|Sari et al. 1998| ) 

p — 2\ nip 



lm^eA'-^\'-^T. (17) 
\p-lj rUe 

The solution to the blast wave dynamics found in the previous Section is valid if the synchrotron 
loss time scale for electrons at the low-energy cutoff 7m is shorter than the dynamical time scale, 
which is comparable to the time elapsed after the explosion. In this case, electrons injected prior 
to any given time will have established a power-law spectrum with index 2 with a low-energy 
cutoff at 

6 7rmeC 



1/2 

(iSari et al. 19981) , where the magnetic field B is parametrized by a constant fraction of the 



equipartition value, 

B = y^32 7r mpC^es n-extT, (19) 

and 7c < 7m for the electrons to be in the fast cooling regime. It is conceivable that the 
magnetic-field equipartition parameter is actually not constant, but increases with time. This 
would not alter Eqs. ^ - ^ (although it would obviously introduce an additional time-dependence 
through ee), but would expand the period of validity of our results because the electrons would 
remain in the fast-cooling regime over a longer fraction of the early blast wave evolution. 

At any given time, the spectrum of electrons behind the shock will then be given by a 
low-energy cutoff at 7c, a power-law with index 2 for 7c < 7 < 7m and a second power-law with 
index p -|- 1 for 7 > 7m. The synchrotron spectrum generated by such an electron distribution 
(neglecting synchrotron-self absorption) consists of a double-broken power-law given by 

[u/ucY/'^ for V < Uc, 

F,=F,{ (u/u,)-^/^ for u,<u< u^, (20) 

(z^m/l^c)"^/^ (i^/j^m)"^/^ for Um 



- 6 - 



where the characteristic synchrotron frequencies z/j = {nieC^/h) (B / Bcr)"fi with Bcr = 4.414-10"^^ G 
and Fc is the flux at frequency Uc- 

Now, normahzing the total number of swept-up, relativistic electrons in the shock to 
Ne = {4:/3)TTnextR^, and the total radiative power to L = (4/3) Tcctt {B^/[8tt]) Ne{-f) -y"^ d'y , 
we find 

„ B Bcr h „ „n 

6b vr dj^ TJieC 

where = 10^^ cm is the luminosity distance to the GRB source. 

We can now use Eqs. |l^ - |2l| to find the relevant scaling laws for the flux normalization Fc 
and the break frequencies t'c and I'm- Denoting by the time in the observer's frame in days, we 
find 

Fc = fF. d^i ef (^^y" To 4" Jy (22) 
fp, = 2.51 • io(-56.o+i5.4 0t+55.3</.B) _ g^)<^t - e)-('^'+<AB), (23) 

'^m = J^A^ f^r" To tf Hz (24) 



with 



(1+z) ^ V 



with 



and 



= 3.74 • 10(12.0+15.4 xt+55.3xi.) (g _ (2 - e)-fe+x^.) , (25) 



-3/2 p52\^- - 



nl- t^' Hz (26) 



with 



fy = 2.40 • 10^^'^-°+^^-^'^*+^^-^'^^) (8 - e)^+'^* (2 - e)"(^+'^*+'^-^\ (27) 

The power-indices of these relations are given in Table |. At a given observing frequency Vobs 
the flux will decay according to power-laws in time, depending on the part of the synchrotron 
spectrum containing the observing frequency at a given time. For < i/c we have Fy oc v^^^ t^" 
with 

u;o = ct>,-- = -^--. (28) 
For Uc < I' < I'm have Fy oc i/^^/^ ^^^i ^^vith 

.,=^„ + ^£ = -2i±i. (29) 
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For u > Um, F,y oc V'Pl'^t'^^ with 

2(l + 6) + 6(p-l) 
a;2 = </)t + y + ^— Xt = . (30) 

Eqs. and ^ are easily inverted to find the sweep-through times and tm of the break 
frequencies Vc and Vm at a given observing frequency, which define the times of breaks in the hght 
curves between the above power-law decay slopes. We have tc oc z/-^/"^* and oc i/-*^/^*. 

From the indices in Tab. || one sees that generally the break frequency decreases more 
rapidly with time than the cooling frequency Uc- If synchrotron self Compton scattering is 
unimportant (see Dermer et al., in preparation), this means that a burst which starts out in 
the fast-cooling regime will become radiatively less efficient, with a time-dependent e, after both 
break frequencies have become equal. The time irf,cm (in days) at which this happens, and which 
terminates the phase of applicability of our results, is given by 

( 33.8 IK/11 I 221 \ 12 8 

td,cm = 10^-— -''■''+14+^^(8 - e)-—^ (2 - e)—^ 

^)-r-*#(^)-, (31, 

where Tcm = (8 -|- 2e)/(8 — e). After this transition, the radiative efficiency will steadily decrease as 
e ~ ee (7m/7c)^"^ QModerski et al. 1999| ) until the blast wave evolution approaches the adiabatic 
limit (e 0) in which, e. g. the expressions given by Sari et al. (1998) may be applied as long as 
the blast wave is relativistic. 



4. The influence of e > 

Fig. |l] shows the GRB afterglow light curves for a standard GRB (see figure caption for 
parameters) at infrared, optical, and soft X-ray frequencies during the relativistic and fast-cooling 



phase of the blast wave (Eqs. ^ - 26 are no longer applicable if the blast wave becomes 
non-relativistic or is in the slow-cooling regime). The figure illustrates that there are very 
pronounced differences between the light curves at a given frequency even for slight changes in 
the radiative efficiency. In particular, the rise/decay slope wqi when the observing frequency u is 
in the Fy oc v'^^'^ part of the spectrum, and the cooling-break sweep-through time tc are seen to 
depend very strongly on the radiative regime. 

The dependence of the temporal rise/decay slopes on the radiative regime, as given by 
eqs. |2^ - 13^ is illustrated in Fig. § for our standard set of GRB parameters. Both the slopes u)q 
and UJ2 (corresponding to the v^f^ and the v~^^'^ parts of the spectrum, respectively) are rapidly 
decreasing with increasing radiative efficiency e. We also see that the typically observed power-law 
decays of X-ray and optical afterglows, F^, oc with 1.1 ;^ a ^ 1.4 (e. g., posta et al. 1997 , 



Piro et al. 199S , in 't Zand et al. 1998 , Galama et al. 1998a , Diercks et al. 19981) , if produced 
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by a blast wave in the fast-cooling regime, imply additional constraints on a rather hard electron 
injection spectrum if one allows for a finite radiative efficiency e ^ 0.1. 

In Fig. ^ we plot the sweep-through times of the break frequencies Vc and Vm for observations 
at optical, infrared, and soft X-ray frequencies for our standard burst parameters as a function 
of the radiative efficiency. In particular for low e, both break times depend very strongly on the 
radiative efficiency. Thus any conclusions drawn from the observed sweep-through times based on 
the adiabatic solution of the blast wave dynamics, need to be taken with caution. Note that in 
calculating these break times as well as the light curves shown in Fig. |l], we have set eg = e, which 
is the reason for the sweep-through time tm approaching for e ^ 0. In the adiabatic limit e — > 
this always leads to tm < tc, implying Um < Vc-, where our solution is not applicable. The figure 
also shows the transition time td^crm at which Vc = Vm- 

Fig. ^ illustrates an example of how the frequency-dependent GRB afterglow light curves are 
expected to change between the early afterglow phase investigated in this paper, and the later, 
most probably adiabatic phase. For the example shown in the figure, we have assumed that the 
blast wave starts out in the fast-cooling regime with e = ee = 0.8. The other parameter values 



are given in the figure caption. According to Eq. 31, our solutions are valid until ~ 35 h into the 
afterglow phase. The following transition phase, in which the blast wave still has a non-negligible, 
but steadily decreasing radiative efficiency, can up to now only be treated numerically. In the late 
afterglow phase, the light curves asymptotically approach the adiabatic evolution, for which we 
are showing the results of Sari et al. (1998). 



5. Summary and conclusions 

We have presented an analytic solution to the hydrodynamic evolution of a relativistic blast 
wave and the broad-band synchrotron spectrum of electrons accelerated behind the shock front in 
general radiative regimes under the assumption that the electrons are in the fast cooling regime 
and the electron acceleration efficiency eg and the magnetic- field equipartition parameter are 
constant. We have shown that in this regime small deviations from a perfectly adiabatic evolution 
will have strong, potentially observable effects on the temporal evolution of the early afterglow 
radiation. In particular, the rise/decay slopes of the monochromatic flux and the times of expected 
breaks in the light curves due to spectral break frequencies sweeping past the observing frequency, 
are significantly affected by a non-zero radiative efficiency. 

Our predictions will be particularly important for the interpretation of rapid optical - X-ray 
follow-up observations of early (~ a few hours after the burst) GRB afterglows, which might 
become possible with much improved data quality with the upcoming HETE II satellite and 
with the planned SWIFT mission and with high-sensitivity, high-spectral-resolution follow-up 
observations by the new generation of X-ray telescopes such as Chandra or XMM. By comparison 
of the decay laws of the afterglow radiation in different spectral regimes and of the shift of the 
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break frequencies with time with our eqs. 28 - ^ and the indices given in Table ffl, the radiative 



regime of the blast wave can be determined unambiguously. Having fixed the radiative efficiency e, 



eqs. and|2^ can be used to determine e^, Fq, and no, if one has an independent estimate of 

the total energy £^2- If such an estimate is not available, one would need an additional equation, 
which can be provided by the location of the synchrotron self absorption frequency, in order to 
solve the system of equations for all four free parameters (e^, Fq, no, and £^52)- 

This work was supported by NASA through grand NAG 5-4055 and the Chandra Postdoctoral 
Fellowship grant number PF 9-10007, awarded by the Chandra X-ray Center, which is operated 
by the Smithsonian Astrophysical Observatory for NASA under contract NAS 8-39073. 
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Table 1: Indices for time, energy, and density scaling of the flux Fc, and the break frequencies 



and fc 



Fc 



= -ih xt = i^t = -2iEf 
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Fig. 1. — GRB afterglow light curves in different radiative regimes at infrared (upper panel), optical 
(middle panel), and soft X-ray (lower panel) frequencies. Parameters: E^2/^ = 100/(47r), no = 100, 
Fo = 300, CB = 10-^ ee = e,p = 2.5, z = 1. 
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Fig. 2. — The temporal decay index LOi in the three different spectral phases as a function of the 
radiative efficiency e. The index L02 (decay in the u'p/'^ part of the synchrotron spectrum) is shown 
for two different values of the electron injection spectral index p. 
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Fig. 3. — The sweep-through times of the coohng and the Vj^ break as a function of the radiative 
efficiency e for three different observing frequencies. Parameters: E^2/^ = 100/(47r), no = 10, 
To = 100, es = 0.2, eg = e, p = 2.5, z = I. Also shown is the time tcm, at which the transition 
from fast coohng to slow cooling occurs (vc = i^c)- The cooling-break sweep-through times tc for 
Ugjjg = 3 ■ 10^^ Hz and 10^^ Hz are of the order ^ 1 s and thus irrelevant since this is shorter than 
the deceleration time scale. 
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Fig. 4. — Transition of the GRB afterglow light curves at infrared, optical, and soft X-ray 
frequencies from the fast-cooling, quasi-radiative regime to the slow-cooling, adiabatic regime. 
Parameters: £^2/^ = 100/(47r), no = 100, Tq = 300, €3 = lO'^, eg = e = 0.8 for t < tcm, 
p = 2.5, z = l. 
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